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Abstract 

Standard derivations of the functional integral in non-equilibrium quantum field 
theory are based on the discrete time representation mna. In this work we der i ve 
the non-equilibrium functional integral for non-interacting bosons and fermions using 
a continuum time approach by accounting for the statistical distribution through the 
boundary conditions and using them to evaluate the Green’s function. 

Non-equilibrium quantum field theory has its foundations on the closed time contour. The 
contour, introduced by Schwinger PISI» is used to evaluate expectation values and replaces 
the real time axis of the equilibrium theory. This contour was developed in statistical quan¬ 
tum mechanics by Kadanoff and Baym [5j and its formulation was made simpler and more 
transparent by Keldysh |6j. Keldyslks formulation is even more transparent in its functional 
integral analogue, developed by Kamenev [7]. This formulation, used also in Refs. mea, 
accounts for the statistical distribution as an off-diagonal component of the inverse Green’s 
function in the discrete time representation. The Green’s function is calculated by matrix 
inversion and the continuum limit is eventually taken. Instead, in this work we adopt a contin¬ 
uum approach: we derive the functional integral by accounting for the statistical distribution 
through the boundary conditions and we evaluate the Green’s function as a solution of a 
differentia! equation with boundary conditions derived from those of the functional integral. 

The Central object of non-equilibrium quantum functional integral is the partition function 
along the closed time contour, 


Tr 


Z = 


U (t h tf)U ( tf,ti)p 


( 1 ) 


Tr [p] 

t j? A 

where p is the initial distribution, lA{tf,ti) = Te -1 ^ dtH ^) (fa — p j s the time-evolution oper- 
ators from the initial time ti to the final time tf and T is the time-order operator [Tj. The oper- 
ators U and U (tj, t/) are referred to as forward and backward time-evolution operators 
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|T]. This partition function lias meaning only if there are different externa! sources on forward 
and backward time-evolution operators. Tliese sources are used to generate correlation func- 
tions and are usually set to zero at the end of a calculation Jl]. When the externa! sources are 
the same on forward and backward time-evolution operators, U ( ti,tf)U {tf,tf) = 1 and the 
partition function is identically Z = 1. However, since the following derivation is independent 
of sources, we omit them, keeping in mind that they should be present. 

We start by constructing the functional integral representation of Z for a non-interacting 
single-level System of energy £ populated by bosons or fermions and we will extend it to many- 
level systems at the end. The Hamiltonian is H{t) = ea^a, where a and a' are bosonic or 
fermionic creation and annihilation operators, satisfying the canonical commutation relation 
[a, a'] + = 1 for bosons or anti-commutation relation [a, a^]_ = 1 for fermions. To derive the 
functional integral representation of the partition function, we consider the coherent States 
of bosons and fermions, (p and (p, respectively associated to a and a t through the relations 
a|0) = (p \(/>) and (cf>\aP = (p((p\, where (p and (p are complex numbers for bosons or Grassmann 
number for fermions. Inserting the resolution of identity of coherent States between the three 
operators in the numerator of Eq. Q and using the coherent state representation of the 
trace, the partition function becomes, 

(4,- \U(U,t f )\4 f ) e-*!*’(4, \U (t f , U) |#> , 

where d[0, <p\ = d(Re0) d(Im0)/7T, ( = +1 for bosons and d[</>, <p\ = d(pd(p, ( = — 1 for 
fermions. We used the superscript + and — for initial states acting respectively on forward 
and backward time-evolution operators. As opposed to the initial states, \(pf) and |</r“), the 
final state, | </>/), acts on botli forward and backward time-evolution operators. 

The partition function assumes a simpler form for initial states defined by the distribution, 


/5 = /W = 


n 


1 + C n 


(3) 


where n is the average occupation of the level, 

n = (a)a) p = 


-/atM -AhM 


Tr 


(4) 


distribution |3 
P = e 


■ , 

T IS t 


Here, n > 0 for bosons and 0 < n < 1 for fermions. An example of a state described by 
is the thermal state of a non-interacting system, where n = (e^r~ — C) -1 ; 
;he Boltzmann factor, T is the temperature and /r is the Chemical potentia!. 
Using distribution (3j) and the identity = e p ^' [T], the integrand in the first line 

of Eq. (2) becomes I {(pt\p\C ( PPP') = and integration over (pf enforces the 

constraint (pf = fpf\ ■ Integrating over (pf the partition function ([2j becomes, 


^ = (!-Cp) C / d[<pi,(pi]d[<pf,<p f ]e ^{(pi\U{ti,t f )\(p f )e ^ f ((p f \U {t f ,U) \(p(pi) , (5) 
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where we remove the label — from the initial state. The pre-factor of the integra! is the 
normalization of the partition function, Tr [p\ = (1 — £p)7 We express the expectation value 
of the time-evolution operator in the f un ction a! integra! form, 

<Ktf)=fy 


Here, the functional integra! measure and the action in the exponent are dehned by splitting 
the time interval tf — ti in N smaller intervals St — (t/ — ti)/N and taking the symbolical 
limit, 

N—l 


v [& 0 ]= 




n= 1 
N 


d t(f>(t) (. id t — e) 4>(t) = lim St 

N—toc ^ J 


n=1 


.7 Vn - Vn-1 T , 

-. ^^Pn^Pn—1 


St 


(7) 


where 0 O = and = (j)f. Integrating over initial and final helds, we arrive at, 


f<j> ( tf)=4> + {tf ) t 

z = (i - c pf / v$ + ,<f>~] ^ 

J <t> + (ti)=cp<t>~ (u) 


where + and — fields are associated to forward and backward time-evolutions and the minus 
sign in front of the second term in the exponent comes from the inversion of the integration 
limits. The integration measure has been redehned to contain the additional integration over 
the initial and final variables, 4>o , 4>o and (f)^, (fy, which are 4>~(ti ), 4>~(ti) and (p + (tf), 4> + (tf ) 
in continuum notation. Tliis partition function is the sum over ali the possible conhgurations 
of + and — fields, related only through the boundary conditions in distant past and future, 
with a weight that depends on the difference of their actions. The stationary conhgurations 
satisfy 0 + (f) = 4>~(t) and 4> + {t) = 4>~(t). These conhgurations are the most important as, 
among the individual conhgurations, they carry most of the weight of the partition functioiQ 
As a consequence, the quantity 0 + (f) — 0 _ (t) measures how far away conhgurations are from 
the most important ones. 


Bosons 

The goal of the remaining derivation is the evaluation of the Green’s function and the pro- 
cedure is different for bosons and fermions. We start from bosons. In analogy to equilib- 
rium quantum held theory |2j, we associate stationary conhgurations, 0 + (f) = 4>~{t ) and 
0 + (f) = (p~{t), to classica! ones. Motivated by tliis identihcation and following Ref. [Tj, we 
dehne a classica! held, 0 cl (f), and a quantum held, 0 q (f), as 

I 0 d (*) = 75 (<^ + W + 0“(*)) > / g v 

\nt) 

1 We note that conhgurations satisfying cf> + = —<f>~ and </> + = —<jr are stationary, as well. However, 
contrary to the previous ones, they may not be stationary for an interacting system [T]. 
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and the same for the complex conjugate. Tliis transformation is known as Keldysli rotation 
In terms of these new fields, classica! configurations have </> q (t) = 0 and, as we will 
see later, </> cl (t) satisfies the classical equations of motion. Moreover, the quantum held 
4> q (t) measures how far away configurations are from classical ones. In these new fields, the 
partition function becomes, 

-i f) = 0 f f _ r 1Q fl 

z = - f dt * a wi G * ] 



where a, f3 = cl, q and we use the conventional sum over repeated indices. The inverse Green’s 
function is, 


[GnT\t) = 


0 

id t - £ 


id t 


0 


( 11 ) 


The inverse Green’s function depends on the initial average occupation, n, because, by con- 
struction, it is an operator whose domain is the set of fields, <p a (t), that satisfy the boundary 
conditions of Eq. © and the initial boundary condition depends on n. Integration over the 
quantum fields, 4> q (t) and 0 q (t), enforces the classical equations of motion, (id t — e)0 cl (t) = 0 
and (—idt — £)</> cl (t) — 0. This fact originally motivated the name “classical field” [lj. 

It is worth to pause for a moment and discuss the initial boundary condition in the 
partition function (10). This condition relates the values of initial classical and quantum 
variables, </> cl (tj) and 0 q (tj), and initial occupation number, n. Taking the modulus, the 
equation reads, 

= (i2) 


Since for bosons n > 0, for a hxed </> cl (tj) the quantum variable is initially bounded, 0 < 
|</> q (fj)| < I0 cl (^i)l (A more restrictive condition applies in the final state, 0 q (t/) = 0). When 
n 1, the level occupation is microscopic and 0 q (tj) ~ 0 cl (fj). Instead, when n 1, the 
level occupation becomes macroscopic and |0 q (tj)| -C |0 cl (t/)|. This, as expected [|2], means 
that the larger the initial occupation of the level, the closer the initial state of the system is 
to a classical one. These results are summarized in Fig. [T[ 

We proceed and calculate the Green’s function, G a/3 (t,t') = —i(0"(f)0 /3 (f / )), wliich satis¬ 
des the differentia! equation, 


f 0 idt -e \ f G clc \t,t') G clq (t, t') \ _ f 5(t-f) 0 \ 

V idt-e 0 J n { G qcl (t,t ') G qq (t,t') J ~ { 0 6(t~ t') J ’ 


where we used the matrix representation (11) of the inverse Green’s function. The general 
solution of this equation is, 


G af3 (t,t') = -ie-' ie{t - t ' ) f 


a (b + 1) 6 (t — t') + bO (t' — t) 

(c + 1) 6 (t — t') + c6 (t' — t) d 


(14) 

The constants a, 6, c and d are hxed by the boundary conditions in (10). For this, we multiply 
each boundary condition of the partition function (10) by 4> a (t'), take the average and use the 
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Figure 1: Example of a configuration of the modulus of quantum and classica! fields, |</> cl (t)| 
and |0 q (£)|. The smaller |</> q (t)| the closer is the conhguration to a classica! one. At the final 
time, t = tf, \4> q (tf)\ = 0. Instead, at t = ti, the initial value of the quantum field, |</> q (tj)|, is 
related to the initial value of the classica! field, |0 cl (fj)|, through Eq. (12). For a fixed value 
of |0 cl (ii)|, \4> q {U)\ —* \4> cl (ti)\ for n —>• 0 and |0 q (t;)| —* 0 for n —* oo. 


dehnition of Green’s function, G a ^(t, t') = —i(0 a (f)0^(f / )), to obtain the boundary conclitions 
of the Green’s function, 

G qa (t f ,t') =0, 

w ’ ’ ’ (15) 

G cla (ti,t') = —(1 + 2n)G qa {ti,t'). 

The first boundary condition fixes the constants of G qa (t, t 1 ) to c = — 1 for a = cl and d = 0 
for a = q, where we use the fact that the first and second theta functions in G qcl (tf,t') 
are respectively equal to 1 and 0 as ti < t’ < tf. Using G qa (ti,t') in the second boundary 
condition, the constants of G cla (t, t ') are fixed to a — 1 + 2n for a = cl and b = 0 for a = q0 
Substituting these values in Eq. 0 , the Green’s function becomes, 




G K (t,f ) G R (t,t') \ 
G A (t. t') 0 ) ’ 


(16) 


where, 

G R (t , t ') = -W (t - t') e _ie ( t_t '), 

G A (t,t') = i9 (t' - t) (17) 

G K (t,t') = -i (1 + 2 n) 


The Green’s functions (|17|) are retarded, advanced and Keldysh Green’s functions, in agree- 

l)" 1 - 


ment with Ref. |TJ for a thermal occupation of the level, n = 


g —M 
T 


Retarded 


2 To fix the constants we did not consider t' = U,tf because this would lead to equal time Green’s functions 
and the theta functions in G qc l (t, t 1 ) and G clq (^,^ , ) have not been regularized at equal times, t = t', yet. 
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and advanced Green’s functions are stili undefined at equal times because the theta func- 
tion, 9(t), needs to be regularized at t — 0. First, using the Fourier transform, G R,A (cu) = 
(cu — Uq + i0 ± ) 1 , we liave, 

G R (t, t ) - G A (t, t) = j ^ [G R (w) - G A (w)] = -i. (18) 

Second, (</> cl (t)</ q (f)) = — ((/> q (t)</> cl (t))' 1 ' which leads to G R (t,t) = [G a (t, i)] 1 ", where f denotes 
complex conjugation and time transposition. It follows that the regularization is 9(t = 0) = 
1 / 2 . 


Fermions 


The derivation is different for fermions, since they do not have a classical analogue. Instead 
of classical and quantum, we use the labeis 1 and 2 through the transformation jTj, 


M*) = ^ O + 0) + 0 (*)) (M*) = ^ (0 + (f) -0 (*)) 

M*) =^(<l> + (t)-(i>-(t)), [Mt) = 75 (0 + W + <i>~(t)) ■ 

In these new variables, the partition function reads, 

where a,b = 1,2 and the inverse Green’s function is, 

id t — e 0 


[g;T‘w = ( 


0 id t — e 


(19) 


( 20 ) 


( 21 ) 


The Green’s function, G ab (t,t') = — i(</> a (t)0 b (t / )), satisfies the differential equation, 


fidt-e 0 \ ( G n (t,t') G 12 (t,t') \ _ f 6{t-t) 0 \ 

V o id t -£ J_\G 21 (t,t') G 22 (t,t') J \ 0 5(t-t’))' [ ’ 

for which the general solution is, 


Qab^. ^ _ _| e -i/ (a + 1 ) 0 (t — t') + ad ( t' — t ) 


(d + 1) 9 (t — t') + d6 ( t' — t ) 


(23) 


We fix the constants a, 6, c and d using the boundary conditions of the Green’s function, 


G 2a (t f ,t') = 0, 

G la (t h t’) = -(l-2n)G 2a (t h t’). 


(24) 
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obtained from the partition function (20) in the same way we did for bosons. The boundary 


conditions at t — tf fix c = 0 and d — — 1 and the ones at t — ti fix a — 0 and b = 1 — 2 n. 


Substituting these values in Eq. (23), the Green’s function becomes, 


G ab {t,t') = 


G K (t,t') G K {t,t') 
0 G A (t,t') 


(25) 


where, 

G R (t,t') = -i 9 (t - t') e-^-^ , 

G A (t, t') = i 0 (t'- t) e-*^, (26) 

G K (t,t') = -i (1-2 n) 

As in the bosonic case, we have retarded, advanced and Keldysh Green’s function, in agree- 
ment with Ref. [1] for a thermal occupation of the level, n = (e - ^ + 1) _1 . The regularization 
of the theta functions are the sanie as for bosons. The only difference between Eqs. and 
(26) is the minus sign in front of n in the Keldysh Green’s function. 


Many levels 

In the previous pages we derived the Green’s function for Systems with a single level. In tliis 
section, we extend the results to systems with many levels. In a System with many levels, 
for every level i, we consider the annihilation and creation operators di and a \, satisfying the 


di 


= 5 , 


ij- 


The single-level non-interacting 
is a Hermitian matrix and 


commutation or anti-commutation relations [a t , a- 
Hamiltonian H = ea)d, is replaced by H = d\eijdj where, e iv 
the average particle number, h, is replaced by the one-body density matrix = (a|dj), 
a Hermitian matrix, as well. For clarity, we employ the bold notation for vectors, a = 
(..., dj_i, di, d* + i, ...), and matrices, e and n. To proceed, we diagonalize the density 
matrix through a unitary transformation, so that the initial distributiori p is the product of 
single-level distributions ([3]), derive the partition function and transform back to the original 
density matrix. Associating the coherent states vectors (p and cp to the operators a and 


a 


t 


the partition function of bosons and fermions are respectively (10) and (20) with bold 


notation and boundary conditions given by, 


(p q {t f ) = 0, 

cp q [ti) = -(/ + 2 n T )- l (p c \ti ), 


and 


02 (tf) = 0, 
(pi(U) = -(T 


2 n T )cp 2 (U ), 


(27) 


(28) 


where I is the identity matrix over the many level indices and T denotes transposition. The 
procedure to calculate the Green’s function is the same as before with the difference that the 
columns on wliich we apply the boundary conditions are over the many-level indices, as well. 
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The Green’s function are identical in form to (16) for bosons and (25) for fermions, but now 
each component is a matrix over the many-level indices given by, 


G R (t,t') = -i 9 {t-t') e-^-^ , 

G A (t , t') = id (t' - t ) e - ie ( t-t '), (29) 

G K (t,t') = -ic ie( ' l ‘ ] (I + 2(n T ) e i£{t '- u) . 


The transposition of the many-body density matrix is due to the fact that the order of creation 
and annihilation operator in the definitions of Green’s function and one-body density matrix 
are exchanged. 


Conclusions 

We saw how to derive the non-equilibrium functional integra! for bosons and fermions by 
accounting for the initia! distribution through the boundary conditions of the functional 
integra! and without resorting to a discrete representation to find the Green’s function. Using 
tliis approach, not only is the derivation shorter and simpler than the discrete approacli but 
highlights the properties of the classica! configurations. Moreover, it is easy to extend the 
procedure to many-level systems with non-diagonal density matrices. 
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